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Abstract 

Starting from the first principles of nonrelativistic QED we have derived the system of 
Maxwell- Schrodinger equations, which can be used for theoretical description of atom op- 
tical phenomena at high densities of atoms and high intensities of the laser radiation. The 
role of multiple atomic transitions between ground and excited states in atom optics has 
been investigated. Nonlinear optical properties of interacting Bose gas are studied: for- 
mula for the refractive index has been derived and the polariton spectrum of a condensate 
interacting with an intense laser field has been investigated. 
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1 Introduction 



In last decades the problem of interaction of the laser radiation with ultracold atomic gases has 
attracted a lot of attention. With the aid of the laser radiation one can manipulate the center- 
of-mass motion of ultracold atoms and one can observe different wave phenomena with ultracold 
atomic beams. After the experimental realization of Bose-Einstein condensation, which allows 
to create rather dense atomic systems, the problem of interaction of photons with ultracold 
atoms has reached a new stage of development. Acting with the laser on a dense atomic sample 
one can induce nonlinearities in the behavior of the matter field caused by dynamical dipole- 
dipole interactions. This provides a possibility to create atomic solitons of different kinds U, 
to change dramatically the effective scattering length of the condensate , to create nonlinear 
beam splitters ||, vortices M, photonic band gaps and defet states in a condensate [Q. 

In recent years different approaches to the description of interaction of ultracold atoms in 
the field of optical radiation have been suggested and different aspects of the phenomenon 
have been considered. The properties of the laser radiation modified by atomic dipole-dipole 
interactions were investigated [§, 0, §[ and it was shown that they can be described by the 
refractive index which is governed in the linear case, when the light intensity is low enough, by 
the Clausius-Mossotti relation known from classical optics if we neglect quantum correlations. 
If the quantum statistical correlations are taken into account, the formula for the refractive 
index contains additional terms defined by a position dependent correlation function || |7], ||. 
Nonlinear optical properties of noninteracting Bose gas were studied as well H . 

The modification of the properties of the laser radiation should have a back influence on the 
behavior of an ultracold atomic ensemble (for example, the motion of atomic beam). An attempt 
to consider this back influence was undertaken by several authors |], |3|, [7], [H| [11], |12[ . In the first 
works on the subject the two-body interactions were modelled by the phenomenological contact 
potential fllQfl . Later on dynamical dipole-dipole interactions were taken into account within 
the framework of the nonrelativistic electrodynamics |], |||, [7|, [ll|, [12|, [13], [14], [15|]. However, in 
papers by W.Zhang and D.Walls and G.Lenz et al.[Q the averaged polarization of ultracold 
atomic ensemble was computed as a function of the incident laser field, whereas it should be a 
function of the macroscopic or the local field which are different from the external laser field 
due to dynamical dipole-dipole interactions. Wallis [|I2J used the correct form of the equations 
for the electromagnetic field. However, his equation for the matter field was not consistent with 



the equation for the electromagnetic field [13]. This problem was considered also by Y.Castin 



and K.M0lmer |TT[ and J.Ruostekoski and J.Javanainen [0]. But the equations which are used 
in those papers are very complicated because they are written down in terms of the local field 
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and dipole-dipole interactions are presented explicitly in the form of the sum over dipole fields. 
This makes the analysis very difficult and the equations, which govern the time evolution of 
the matter fields, remained unsolved. 
In the papers f 



14] , [jjl a self-consistent quantum theory of atom optical processes has been 



developed. Making use of the Lorentz-Lorenz relation, which allows to simplify the analysis, 
it was obtained the general system of Maxwell- Schrodinger equations for atomic creation and 
annihilation operators and the propagation equation for the laser field which can be used for 
the description of linear and nonlinear phenomena in atom optics of single- and multi-species 
condensates at high densities of the atomic system. However, the treatment in [13[ |14], [15[] was 
restricted by low light intensities. 



In the present paper we shall continue the investigations started in |13|, 14]. Having in mind 
mainly atom optical applications, we shall derive the system of Maxwell-Schrodinger equations 
for the case of high density of the atoms and high light intensity. Nonlinear optical properties 
of the interacting Bose gas will be also discussed. 



2 Heisenberg equations of motion for the atomic opera- 
tors 

We consider a system of ultracold two-level atoms with masses m, transition frequencies u a , 
and transition dipole moments d. We shall describe such a system in terms of matter field 
operators. Let \g) and |e) are the vectors of the ground and excited states of the quantized 
atomic fields. Then the corresponding annihilation operators of the atoms in these internal 
states are (fi g and <p e . Matter- field operators are assumed to satisfy the bosonic equal time 
commutation relations. 

The Heisenberg equations of motion for the atomic operators are derived from the Hamil- 
tonian of the second quantized atomic field interacting with the photons. In the reference 
frame rotating with the frequency ujl = cki of the incident laser field Ej n (r, t), which is as- 
sumed to be monochromatic, and in making use of the electric dipole approximation and the 
rotating-wave approximation we obtain the following dynamical equations for the matter-field 
operators [|13|, [L4| 

ih-Qf = H cm (f) g + H ge (f> e , (1) 
itS^ = H C m4>e - ft (A + ry/2) e + H eg 4>g , (2) 
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where H cm = —Ti 2 V 2 / (2m), A = uj l — cu a — 5 is the detuning of the frequency of the laser wave 
from the frequency of the atomic transition, 8 and 7 are the Lamb shift and the spontaneous 
emission rate of a single atom in free space, respectively. Here we have introduced the operators 
H eg = — dE^ c and H ge = — dEf oc which are responsible for the transitions \g) — > |e) and 
|e) — > \g), respectively. H eg and H ge are related to the operator of the local electric field 
E^ c (r, t). The positive- frequency part of this operator has the form 



EL( r > *) = E m(r) + i £ \/^7^e A c kA (0) exp [ikr -i(u k - uj l ) t] 



+ ( drV x V x P + ( f/ ^ R / c ) e ik L R ; (3) 
J R 

where Vx refers to the point r and the operator CkA(O) corresponds to the free-space photon 
field (vacuum fluctuations). The polarization operator P + = d0^0 e . Note that in equation 
(|3|) a small volume around the observation point r is excluded from the integration. The last 
term in equation (|3|) describes dipole-dipole interaction. Due to the structure of the local field 
given by equation ([|) one can distinguish two kinds of the photons: primary laser photons and 
secondary photons, re-radiated by the atoms. 

Usually in atom optics of BEC one deals with equations for the ground state matter-field 
operator S . Therefore, one has to eliminate the excited state matter-field operator e from 
the system of equations (|1|), (0). In the next section we shall develop a general procedure of 
the elimination of the excited state which is in contrast to previous works (see, for instance, 
H) valid for high densities of the atoms and high intensities of the laser radiation. 

3 Elimination of the excited state 

Let's assume that initially there are no atoms in the excited state. Then we can rewrite equation 
(|2|) in the form 

kit) = -i [\i^- A y H^e-^'Mt'W , (4) 

J 

where A = A + Making use of the identity 

ax 00 / -i\k ak 

e™F(x)dx = — £ L-jL—Fix) , (5) 
a 1~ n a" ox* 

k=0 

which is valid for a no n- vanishing a, we get 



-iZfO-t OO (_-\k f)k r . 

e(r , t ) = V e i3 f**H eg (r, t)MT, t) 

Y v ' ; hA ^ A k dt k [ gK ' ' ; 



(6) 



4 



If H eg does not depend on time and the center-of-mass motion does not give a large contribution, 
we have 

(7) 



H (r) 00 (-i) k B k 



hA t^o A fc dtk 

From the system of equations ([I]),(|2]), with the center-of-mass motion neglected one can derive 
the relation 



dH X (r t) = _ z A(r)i/ eg (r) ~ Hl!^ (rt) n=12 
'V*-'J- h 2 A 2 ~ n ^_ A k dt k ^ 9{ - , h 



dt 



Substituting this relation iteratively into equation (0), we obtain 



hA 



m=0 



# 9 e(r)ff eg (r) 



9 M) 



(9) 



«o = 1, 

1 fci + lfe 2 + l fcm-i+l 

wi = E I! E ■" E 

fc 1= fc 2 =0 fc 3 =0 fc m =0 



, (2m + 1)! 
m\(m + 2)! 



m = 0, 1, 



(10) 



The zeroth order term (m = 0) in equation ([]) corresponds to the transition \g) — > |e). The 
next term (m = 1) corresponds to the transition \g) — > |e) — > |g) — > |e) and so on. Therefore, 
different terms in equation @ describe multiple transitions of the atoms between ground and 
excited states, caused by the influence of the photons. In a typical atom optical situation when 
the electromagnetic field has a form of a standing wave these terms correspond to the processes 
with the momentum transfer from the laser beam to the atoms ±2h~K(m + 1), where K is a 
wave vector of the laser wave in a medium. 

We assume that e = H ge H eg /(h 2 A 2 ) acts only on states for which the series in equation 
(|9|) converges, i.e., \ip) can be decomposed into eigenstates of e whose eigenvalues satisfy the 
condition |e n | < 1/4. Then the result of the summation is given by 



0e(*) 



H eq VT+Te - 1 



HA 2e 

and the equation for the ground state (fl]) takes the form 



fill 



ih^- 
dt 



H, 



hA 



rinY g + — (Vl +46 - l) <Pg(t) 



(12) 

Although we do not see explicitly on the r.h.s. of equations (JTTJ) , ( JT2] ) e , it is still there, 
because E^ c depends on <p e and this dependence, which is given by fl3|), is very complicated. 
Formally we could iteratively substitute equations (|J) and (]TT|) into (0). However, by doing 
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such a procedure we would get an equation for the ground state wave function in a form which 
would be impossible to use. Therefore, the algorithm described in this section does really 
eliminate the excited state only in the case when the dipole-dipole interactions are negligible, 
i.e., when E^ c m Egj. In order to eliminate the excited state in the case when the dipole- 
dipole interactions play an important role it is useful to employ in addition the Lorentz-Lorenz 
relation, which will be briefly discussed in the next section. 

Because we are mainly interested in atom optical problems and want to study the coher- 
ent evolution of the center-of-mass motion of the gas, we shall neglect spontaneous emission. 
This is valid for situations where the absolute values of the detunings are much bigger than 
the spontaneous emission rates and Rabi frequencies |A| 3> 7 , dE^ c /h . In order to do 
this approximation self-consistently we drop in the following the vacuum fluctuations and the 
spontaneous emission rates 7 from our equations. Then we may replace all the operators by 
macroscopic functions. 



4 Local-field correction 

As it was mentioned in the previous section, the solution of equations (^), flTTp, ( |T!?D is a rather 
complicated mathematical problem because these equations contain explicitly the dipole-dipole 
interactions. In many particular situations such a detailed microscopic description of matter 
is not necessary and it is more convenient to consider optical properties of the medium on a 
macroscopic level. This can be done by introducing the macroscopic field E mac (r,t), instead of 
the local field E; oc (r,t) in the equations for the matter fields. 

As in Ref. |16] we can introduce the macroscopic field by imposing the requirement that it is 
a solution of the macroscopic Maxwell equations for a charge-free and current-free polarization 
medium. Using Maxwell equations and the definition of the local field (|3]), we get the following 
relation 

ELM) = E± ac (r,*) + y P ± (r,t) . (13) 

This equation is often called in the literature the Lorentz-Lorenz relation. It constitutes the 
basis of the local-field effects in classical, quantum and nonlinear optics. 



5 Nonlinear matter equation 

Let's consider first a special case of the low light intensity. In this case we can keep only the 
first term (m = 0) in equation ([|), which is linear with respect to the electromagnetic field 
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strength: 

From equations (|T3|), ( |T4] ) we obtain 

0+(r) 



dE+ c (r) 



e (r,t) 



2A 



n=0 



4tt 2 >" 



(14) 



(15) 



where the position dependent Rabi frequency fi + (r) = 2dE+ ac (r)/7i is related to the macro- 
scopic electric field, a = —d 2 /hA is the atomic polarizability. The series in (|T5|) converges, 

'g I 



provided that ^-\a\ \(pg\ < 1. In this case we get 



e (r,t) 



2A 



: -a\<P g (r,t)\' 



(16) 



This adiabatic solution has been obtained in Ref. [13] using slightly different technique. Note 

2 

= 1. However, as it 
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a 



that a singularity occurs in equation (|TjJ) under the condition 
follows from our present derivation, we never encounter this singularity, because in the region 
| 2 > 1 equation (|T6"D is not valid. This was not clear from the derivation given in 
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a 



Ref. 13 



Then substituting ( Jl6| ) in (JJ), we obtain as the result an equation for the ground state 
1131,111 



matter field 



ift. 



dt 



2V72 



2m 



h\tt + (r) 



4A 



l-fa|0 s (r,t)r 



5 ( r ,*) 



(17) 



Now we substitute the Lorentz-Lorenz relation ( p!3[ ) into equations (pTT|) , (|12|) and keep the 
terms up to the order 1/A 3 . As a result we get the following equation for the ground-state 
wave function 



ih- 



Ot 



H, 



cmfg 



+ 



h\n- 



i2 r 



4A 



8tt 



1 H a \4>a\ +3 — a 



An 



]n+| 

4A 2 



21 



J 9 • 



Varying the density of atoms \<§> g \ , the light intensity, which is proportional to , and 
the magnitude and the sign of the detuning A, one can change the effective potential for the 
ground-state matter field. 

It is interesting to compare the nonlinear terms in equation ( |T8D with the terms in the 
equation for the condensate wave function, which takes into account the effects of quantum 



fluctuations [|17]] 



ih- 



dt 



H, 



crafg 



+ 



Anh 2 a 



m 



1 32 fl3 , j. 

1 + 



J 9 ' 



(19) 
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where a is a scattering length. The leading nonlinear terms in equations (0), (0) are defined 
by the quantities 

respectively. The typical orders of magnitude of the parameters for alkali metal atoms are (in 
CGS system of units) d ~ 10~ 18 esu, \a\ ~ 10~ 7 cm, m ~ 10~ 23 g, and we get the estimates 

We\ ~ ^r^lO" 36 erg cm" 3 , \U C \ ~ 10 -37 erg cm" 3 . 
A^ 

Therefore, we see that at |fi+r /A 2 ~ 0.1 U & and U c are of the same order of magnitude. 
Higher order nonlinear terms in equations flUf) and ( |i~9"|) are defined by 

respectively. For the densities |0 9 | 2 ~ 10 14 cm -3 and detunings |A| ~ 10 8 Hz, we have 

IfJ+l 2 

\ V e\ ~ 'J 10 38 erg cm 3 , \V C \ ~ 10 39 erg cm 3 . 
A 2 

Thus, we see that the nonlinear terms in equations (|18D and (|T9|) have the same order of 
magnitude. 



6 Optical properties of the ultracold gas 

6.1 Refractive index 

We substitute equation ([11]) into the definition of the polarization field. This gives us a general 
nonlinear relation between the polarization and the local field: 



>±_„,vT+4e-1 2 . 



a- 



2e 



\UKc- (20) 



In the special case of the low light intensity we can use the adiabatic solution (jig ). Then 
the expression for the polarization takes the form 

P ± = a\<j )g \ 2 Et = xKac, (21) 
where dielectric susceptibility is given by 



a |0, 



|2 



X — An I , ,2 • (22) 
1 - T a\<j) g \ 
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This expression for the dielectric susceptibility x leads to the Clausius-Mossotti formula for 
the refractive index. As it follows from our physical interpretation of the expansion (|9|), the 
Clausius-Mossotti formula corresponds to the quantum transition of the type \g) — > |e). There- 
fore, it takes into account pair interactions between the atoms when one atom emits a photon, 
then this photon is absorbed by another atom and so on. 

If we keep the terms up to the order 1/A 3 , the expression for the polarization takes the 
form 



a \(j) g 



1 



dE 



loc 



Ti IS? 



± 

mac i 



where 



X = a\(j> 9 



An 2 
1 + — oc\(j)g\ + 



4?r u I 2 
y« \<Pg\ 



|dE+ I 

I mac I 



(23) 



(24) 



3 "'" yi ' V 3 "'" yi J h 2 A 2 
Dielectric susceptibility is a rather important parameter, because it describes the propa- 
gation of the laser radiation inside a medium. In most of the practical situations the electro- 
magnetic processes are much faster than the center-of-mass motion of the atoms. Therefore, \ 
can be considered as a time-independent quantity. Let us assume in addition that the spatial 
variations of the atomic density are not very large, such that Vx — ► 0. Then divE^ ac rs 0, and 
we have the following Helmholtz equation for the macroscopic electric field 



V 2 E± ac + kjn'E] 



2-™± 
L lh ^ 'mac 







(25) 



with the refractive index n given by 



n 



1 + 4irx = 1 + 27ra 



7T 
— ( 

3 



5 /47T 

+ , 

8 V 3 



l + -a|0 fl r+ o l-5-«iv ff i 



IdE: 



i2 -| 



h 2 & 



(26) 



which corresponds to the Kerr-type optical nonlinearity. 

Equations (|TE|), (P5|), ( ^B] ) can be considered as an atom optical analogue of the system of 
Maxwell-Schrodinger equations used in quantum and nonlinear optics. In general they have to 
be solved in a self-consistent way. 



6.2 Polariton band gap 

In this section we investigate the polariton spectrum of a condensate interacting with an intense 
laser field. Let's consider a special case of dilute atomic gas. In this case one can neglect the 
local-field correction in equation (p~3|) and we have E Zoc (r, t) ~ E^ ac (r, t). In order to simplify 
the analytical analysis we assume that the intensity of the laser radiation is not too high. Then 



9 



from Maxwell equations and equation 
index 

2t \ -i R(uJ — UJ 



n\uo) = 1 



we get the following expression for the refractive 
4vrd 2 \d) n \ 2 



R 



(cu -u o y + /4 n 
which contains saturation effects, owing to the intense electric field. 
The polariton spectrum uj(K) can be obtained from the equation 

UJ 



(27) 



K 



-n(ui) 



(28) 



As it follows from equation fl2"?|) n{oS) is imaginary for the frequencies uj within the range 

r-Jr 2 - \Q+\ 2 



uj + 



< UJ < UJq + 



R+\ R 2 - 10+1 



R > 



(29) 



In this case equation (P8|) has no real solutions. This means that light can not propagate in 
a condensate, i.e., we get a polariton band gap with the width Aq = J R 2 — |^ + | 2 - The gap 



width Aq increases with the density of the atoms in the ground state \(j> g \ and decreases with 
the laser intensity. In the case of the low light intensity (|fi + | <C R) we get a well-known result 
Aq = R |TB, If R < |^ + |, equation fl2"8p has only real solutions. This means that in this 
case the polariton gap disappears and the medium is always transparent. 

Note that our analysis does not take into accout effects caused by the influence of the 
spontaneous emission. As it was discussed in Ref. [Tj|, the spontaneous emission smears out 
the polariton spectrum, providing accessible polariton states inside the gap. 



7 Conclusion 

Starting from the microscopic model and making use of the multipolar formulation of QED, 
we have derived the general system of Maxwell- Schrodinger equations for atomic creation and 
annihilation operators and the propagation equation for the laser field. It describes the modifi- 
cation of the properties of the external off-resonant laser radiation in a medium due to dipole- 
dipole interactions and the influence of this modification on the center-of-mass motion of the 
ultracold atoms as a single dynamical process. The system can be used, for instance, for the 
self-consistent analysis of linear and nonlinear phenomena in atom optics at high densities of 
the atomic system and high intensities of the laser radiation. 

A general procedure of the elimination of the excited state has been developed. The anni- 
hilation and creation operators of the excited state for large detuning are represented in the 
form of a series expansion in powers of the inverse detuning, which corresponds to multiple 
transitions between the ground and excited electronic states of the atoms. 
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Optical properties of an interacting ultracold Bose gas are studied: formula for the intensity- 
dependent refractive index is derived and the polariton spectrum of a condensate interacting 
with an intense laser field is investigated. 

Acknowledgments 

This work has been supported by the Deutsche Forschungsgemeinschaft and the Optikzentrum 
Konstanz. One of us (K.V.K.) would like to thank also the Alexander- von-Humboldt Stiftung 
for financial support. Helpful discussion with G.V.Shlyapnikov is gratefully acknowledged. 



11 



References 

[1] Lenz, G., Meystre, P., and Wright, E., 1994, Phys.Rev.A, 50, 1681. 

[2] Fedichev, P.O., Kagan, Yu., Shlyapnikov, G.V., and Walraven, J.T.M., 1996, 
Phys. Rev. Lett, 77, 2913; Bohn, J.L. and Julienne, PS., 1997, Phys.Rev.A, 56, 1486; 
Cornish, S.L., Claussen, N.R., Roberts, J.L., Cornell, E.A., and Wieman, C.E., 2000, 
Phys. Rev. Lett, 85, 1795. 

[3] Zhang, W., and Walls, D., 1994, Phys.Rev.A, 49, 3799. 

[4] Marzlin, K.-P. and Zhang, W., 1998, Phys.Rev.A, 57, 3801; Marzlin, K.-P. and Zhang, 
W., 1998, Phys.Rev.A, 57, 4761; Dobrek, L., Gajda, M., Lewenstein, M., Sengstock, K., 
Birkl, G., and Ertmer, W., 1999, Phys.Rev.A, 60, R3381; Marzlin, K.-P, Zhang, W., and 
Sanders, B.C., 2000, Phys.Rev.A, 62, 013602. 

[5] Marzlin, K.-P. and Zhang, W., 1999, Phys.Rev.A, 59, 2982. 

[6] Morice, O., Castin, Y., and Dalibard, J., 1995, Phys.Rev.A, 51, 3896. 

[7] Ruostekoski, J., and Javanainen, J., 1997, Phys.Rev.A, 55, 513. 

[8] Ruostekoski, J., and Javanainen, J., 1999, Phys. Rev. Lett., 82, 4741. 

[9] Kasevich, M.A. and Harris, S.E., 1996, Optics Letters, 21, 677. 

[10] Lenz, G., Meystre, P., and Wright, E., 1993, Phys. Rev. Lett, 71, 3271. 

[11] Castin, Y., and M0lmer, K., 1995, Phys.Rev.A, 51, R3426. 

[12] Wallis, H., 1997, Phys.Rev.A, 56, 2060. 

[13] Krutitsky, K.V., Burgbacher, F., and Audretsch, J., 1999, Phys.Rev.A, 59, 1517. 

[14] Krutitsky, K.V., Burgbacher, F., and Audretsch, J., 2000, Laser Physics, 10, 15. 

[15] Krutitsky, K.V., Marzlin, K.-P., and Audretsch, J., 2001, Laser Physics, 11, 449. 

[16] Ghiner, A.V., and Surdutovich, G.I., 1994, Phys.Rev.A, 49, 1313; 1997, Phys.Rev.E, 56, 
6123. 

[17] Braaten, E. and Nieto, A., 1997, Phys. Rev. B, 56, 14745; Anderson, J.O. and Braaten, E., 
1999, Phys.Rev.A, 60, 2330. 



12 



[18] Svistunov, B.V. and Shlyapnikov, G.V., 1990, Zh.Eksp.Teor.Fiz., 98, 129. 
[19] Politzer, H.D., 1991, Phys.Rev.A, 43, 6444. 



13 



